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ABSTRACT 



High spatial resolution low-/ 12 CO observations have shown that the wide-angle outflow seen in the Orion BN/KL region correlates 
with the famous IL fingers. Recently, high-resolution large-scale mappings of mid- and higher-/ CO emissions have been reported 
toward the Orion molecular cloud 1 core region using the APEX telescope. Therefore, it is of interest to investigate this outflow in the 
higher-/ l2 CO emission, which is likely excited by shocks. The observations were carried out using the dual-color heterodyne array 
CHAMP + on the APEX telescope. The images of the Orion BN/KL region were obtained in the 12 CO 7=6-5 and 7=7-6 transitions 
with angular resolutions of 8'.'6 and 7'.'4, respectively. The results show a good agreement between our higher-/ 12 CO emission and 
SMA low-/ 12 CO data, which indicates that this wide-angle outflow in Orion BN/KL is likely the result of an explosive event that 
is related to the runaway objects from a dynamically decayed multiple system. From our observations, we estimate that the kinetic 
energy of this explosive outflow is about 1 - 2 x 10 47 erg. In addition, a scenario has been proposed where part of the outflow is 
decelerated and absorbed in the cloud to explain the lack of CO bullets in the southern part of BN/KL, which in turn induces the 
methanol masers seen in this region. 
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1. Introduction 



radio source /, and/or source n, which may be runaway ob 



Orion BN/KL dBecklin & NeugebaueR 1967; 

iKleinmann & LowL 1 19671) . the most luminous source within 
the Orion molecular cloud 1 (OMC-1), has an infrared (IR) 
luminosity of ~ 10 5 L and has been intensively observed 
in numerous molecu lar lines and also at various wavelengths 
dO'Dell et all 120081) . Orion BN/KL harbors the eponymous 
mole cular "hot core" (H C) with a rich and complex chemistry 
(e.g., iBlake et ail Il987l) . In addition, shock and outflow phe- 
nomena are also widely present in the Orion BN/KL region. 
Near-infrared (NIR) obser vations showed H2 outflows in the 
form of bullets or fingers dTavlor et all 1 1984 lAllen & BurtonL 
fl993h at a wide angle along a northwest-southeast axis. In 
radio observations, the powerful ou tflow was also seen in 
the high-velocity 12 CO line wings dKwan & Scovi M U976I) 
togeth er with some hydroxy 1 (OH) mas ers that trace the shock 
fronts dNorrisi 11984 ICohen et all |2006). Another outflow with 
lower v elocities was detected through proper m otions of H2O 
masers dGenzel et all [l98lt iGaume et all 119981) accompanied 
with SiO masers near t he IR source IRc2 dGreenhill et allll998t 
iDoeleman et afl IT999). and is closely associated with source / 
( Menten & ReidL 1 1 9951; iReid et all 120071) . 

The origin of the wide-angle high-velocity outflow and 
northeast-southwest low-velocity outflow has been frequently 
discus sed, and new debat es have begun with updated observa- 
tions. IZapata et alJ d2009l) reported explosive 12 CO bullets in 
high spatial resolution associated with the H2 fingers, and in- 
dicated that this outflow is caused by disintegration of a mas- 
sive stellar system that happened about 500 years ago. Evidence 
of a close dynamical interaction of a young stellar system is 
also supported by the proper motions of the Orion BN source, 



jects f rom a dynam i cally decayed multiple sys tem (Nissen et al 



2005 



20 lj iGoddi et all l20!Tl iBallv et all I20T1I iRodrfguez et al 
I et ali r 



Gom ez et al ., 2005). Such an event could also produce the 



energy needed for driving the large-scale outflows in the Orion 
BN/KL region (Ball y & Zinnecketl 120051) and may indeed be re- 
sponsible for most of the luminosity in the region because unlike 
other hot molecular cores, the Orion HC (mm/submm continuum 
emission peak) does not appear to contain an embedded mas- 
sive protostellar object c apable of producing the necessary en- 
ergy (Za pata et alll201 H) . In this paper, we present the results of 
the new 12 CO 7=6-5 and 7= 7-6 observations toward the Orion 
BN/KL explosive outflow. 

2. Observations 

The observations were carried out in 2007 November and 2008 
July with the 12 meter APEX telescope on Llano de Chajnantor 
in Chile dGtisten et all [200613, an d were part of the large-scale 
on-the-fly (OTF) mapping in the OMC- 1 co re region (for more 
observational details, see iPeng et all 120121) . The observations 
of the 12 CO 7=6-5 (691.5 GHz) and 7=7-6 (806.7 GHz) lines 
were performed with the dual-color he t erody ne array CHAMP + 
dKasemann et all l2006t iGusten et all 120081) . operating 2x7 
elements simultaneously in the 450 and 350 /urn atmospheric 
windows. The MPIfR Array Correlator Spectromete r (MACS) 
and the fast Fourier transform spectrometers (FFTS. lKlein et al 



1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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Fig. 1. a) Redshifted 12 CO 7=6-5 emission (left panel) at the 15 km s _1 channel (A V = 1 km s _1 ) and the blueshifted 12 CO emissio n 
(right panel) at the 3 km s -1 channel (AV = 1 km s -1 ) overlaid with the SMA 12 CO 7=2-1 bullets taken from Zap ata et al.1 J2009). 
The 12 CO 7=6-5 emission is plotted in dark red and dark blue contours running from 8 to 48 K km s in steps of 4 k km s (1 
cr), and the contours with intensities higher than 48 K km s _1 were skipped for clarity. The redshifted (from +40 to +120 km s ) 
and blueshifted (from -30 to -110 km s" 1 ) 12 CO 7=2-1 emission peaks at different velocities are plotted as red and blue dots. 
The squares represent the positions of the outflow explosion center and the H2 emission peak 1 and 2. The triangles represent the 
positions of the Orion HC (south) and BN (north). The stars mark the positions of the five Trapezium stars (0 l Ori A, B, C, D, and 
E). The positions of the H2 emission peak 1 and 2 are marked as squares. Note that the emission in the southwest of the Trapezium 
stars comes from the OMC-1S region, b) 12 CO 7=6-5 (black) and 7=7-6 (red) spectra toward Orion BN (upper panel) and the HC 
(lower panel). Dashed lines indicate the integrated velocity ranges shown in Figs. [2] and [3] 



2006) were used as backends. Spectral resolutions were ~ 0.42 
km s 1 for the MACS and ~ 0.05 km s" 1 for the FFTS. The an- 
gular resolutions of CHAMP+ are 7'.'4 at 806.7 GHz and 876 at 
691.5 GHz. The maps were centered on the Orion BN source 
(a, 5/2000 = 05 h 35 m 14 s .16, -05°22'21'.'5). During the observa- 
tions, position switching to an emission-free reference position 
at offsets (-500", 0") was used. Mira (o Ceti) was used as a line- 
pointing source in the 12 CO 7=6-5 line, and the overall pointing 
accuracy is estimated to be < 4". All spectra were converted 
to the main beam brightness temperature unit, 7mb = T^/j]mb 
(?7mb is 0.47 and 0.45 for I2 CO 7=6-5 and 7=7-6, respectively), 
using a forward efficiency of 0.95 and beam coupling efficien- 
cies as determined toward Jupiter. All data were reduced using 
the standard procedures in the GILDAS packageQ. 

3. Results and discussion 

In Figure Q] a, the faint high-velocity I2 CO 7=6-5 emission 
around the BN/KL core region shows a clear footprint from 
the finger-like structure revealed by the SMA 12 CO 7=2-1 bul- 
lets (IZapata et all H009), especially in the northwest (redshifted 
components) and north (blueshifted components). In addition, 
Figure [2] shows the very high velocity (up to +100 km s -1 ) 
components in the 12 CO 7=7-6 emission, and the 12 CO out- 
flows are clearly associated with the H2 emission peak 1 and 
2, th e two strongest peaks of vibrationally excited H2 emis- 
sion (ISugai et allll994fclBeckwith et all[l978h . The 12 CO 7=7- 
6 and 7=6-5 spectra toward the Orion HC and BN are shown 
in Figure Q] b, where the extended line wings up to ±100 
km s" 1 are clearly visible. This very high velocity outflow has 
been observed in different 12 CO transitions before (see, e.g. , 
IWilson et all 1200 lURodriguez-Franco et allll999USchulz et all 



2 http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 2. Image of the very high velocity emission in the CO 
7=7-6 transition (smoothed by a Gaussian profile with a width 
of 16'.'4). Red contours are plotted from 30 to 90 K km s in 
steps of 10 K km s _1 (2 cr level), and blue contours are plotted 
from 30 to 130 K km s in steps of 10 K km s (2 cr level). 
The positions of source /, the HC, SMA-1, and source n are rep- 
resented as triangles, and their relative positions are indicated 
in the bottom-right panel. The squares represent the positions of 
the outflow explosion center (EC) and the H2 emission peak 1 
and 2. 



119951: ISchmid-Burgketall 1 19891) . and the morphology of the 
high-velo city 12 CO 7=7-6 emission is similar to the result ob- 
tained by Rodriguez-Fra nco et al.l (1 19991) . where the high-speed 
blueshifted and redshifted bullets seem to be associated with 
and preceded by the H2 emission peak 1 and 2, respectively. 
However, the up per limit of the H? outflow radial velocity is 
~ 34 - 40 km s _1 dNissen et all 120071) . which is only half the ve- 
locity of the high-speed 12 CO bullets. Another noteworthy simi- 
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Fig. 3. a) Very high velocity outflow emission in the 12 CO 7=6-5 transition overlaid with the SMA 12 CO 7=2-1 emission (# syn = 
3'.'3 x 3'.'1, P.A.=-14°. IZapata et all l2009h and the H 2 column density dPeng et al.Ll2012h . The left panel shows the redshifted 12 CO 
emission integrated from +70 to +100 km s _1 (Vlsr), and the middle panel shows the blueshifted 12 CO emission integrated from 
-100 to -70 km s . The 12 CO 7=6-5 emission is plotted in orange and light blue contours running from 10% to 90% in steps of 
10% of the peak intensity. The 12 CO 7=2-1 emission is plotted in thick r ed and blue conto urs from 30% to 90% in steps of 20% of 
the peak intensity. The H2 column density (derived from CO isotopologs. lPeng et al.Ll2012h is plotted in gray contours running from 
10% to 90% in steps of 10% of the peak density (5.8 x 10 23 cm -2 ). The right panel shows the superposition of the high-velocity 
12 CO emission from the left and middle panels, and the beam sizes of APEX and SMA are indicated in the first panel. The squares 
represent the positions of the outflow explosion center (EC) and the H2 emission peak 1 and 2. The triangles represent the positions 
of the Orion HC (south) and BN (north). The stars mark the positions of the five Trapezium stars {8 l Ori A, B, C, D, and E). b) The 
corresponding mean intensities (redshifted and blueshifted) are measured in different position angles (in steps of 45°). Note that the 
SMA 12 CO 7=2-1 intensities were multiplied by 3. 



Table 1. Orion BN/KL explosive outflow parameters 



Component" 


f 


A APEX " 


Asma 6 




Energy'" 




(%) 


(10- 


3 2\ 
pc-) 


(M ) 


(10 46 erg) 


[-100,-70] 


-22 


6.5 


1.4 


0.20-0.24 


1.8-2.2 


[-70, -30] 


-25 


7.7 


1.9 


1.71-2.09 


6.1-7.5 


[+30, +70] 


-26 


8.9 


2.3 


2.61-3.18 


4.1-5.1 


[+70, +100] 


-14 


8.8 


1.2 


0.41-0.50 


2.3-2.8 



Notes. (a) Integrated velocity ranges in Vlsr- {b) The filling factors were 
estimated by comparing the 12 CO 7=6-5 emission sizes (A apex) with 
those of the SMA 12 CO 7=2-1 emission (A S ma) at the contours of 20% 
of the peak intensities. (c) The outflow energy (5-M ol]lflow V 2 ) was calcu- 
lated by using the median velocities for the redshifted and blueshifted 
components. 



larity between the H2 fingers and 12 CO outflows is that the wide- 
angle distribution in the northwestern part is presented both in 
the H2 observations and our data. 

The very high velocity 12 CO 7=2-1 emission (Fig. 0a) also 
shows clear emission counterparts in the APEX 12 CO 7=6-5 
emission, where the redshifted emission lies mostly in the NW- 
SE direction and the blueshifted one lies mostly in the N-S direc- 
tion. The average intensities over different position angles (Fig. 
[3]b) also indicate the similarity between SMA 12 CO 7=2-1 and 
APEX 7=6-5 lines. Interestingly, the dense ridge peaking at the 
Orion HC is located to the east of the explosion center. 

The column density of the outflow can be estimated using 



V( 12 CO) 



3kQ n 



8tt 3 v5jU 2 / 



TmbcIV, 



(1) 



assuming optically thin 12 CO 7=6-5 emission at high veloc- 
ities. The optical depth of 12 CO 7=6-5 decreases to about 
0.1 at 20 km s _1 , estimated from the C 18 (6-5)/ I2 CO (6-5) 



ratios dPeng et all I20T2I) with a r 12 CO]/rC 1 8 01 r atio of 490 
dBoreiko & BetzL 119961: IWilson & Matteuccjl Il992l) . Q mt is the 



rotational partition function, fi the dipole moment of 12 CO, S the 
line strength, and 7 ex the excitation temperature. The upper level 
energy E up /k for 12 CO 7=6-5 is 1 16 K. The beam filling factor 
/ was estimated by comparing the outflow sizes (at the contours 
with 20% of peak intensities) of the 12 CO 7=2-1 interferome- 
try data with the APEX 12 CO 7=6-5 emission sizes, given that 
these bullets were resolved by the SMA beam (3'.' 28 x 3'.' 12). 
The excitation temperature 7 ex of about 60 K was estimated 
from the ratio between 7=2-1 and 7=6-5 in the LTE condition. 
The outflow mass and kinetic energy then can be calculated by 
ad opting a r 12 CQ1/ r H 2 l abund ance ratio of 8 x 10" 5 in OMC- 
1 (IWilson & MatteuccI 119921) . The derived parameters of this 
outflow are listed in Table Q] The total kinetic energy associated 
with this explosive outflow in Orion BN/KL is about 1 - 2x 10 47 
erg, close to the value derived bv lKwan & Scovillei d 19761) with 
a revised distance to the Orion Nebula (414+7 pc, Menten et al., 
2007). The outflow gas with low velocities (-30 < Vlsr < 30 
km s~') is likely contaminated by ambient gas and was ex- 
cluded from the calculation (the filling factor is also difficult 
to estimate due to the filtering effect of the SMA data). It is 
likely that the low-velocity gas only contributes a small frac- 
tion of the total outflow mass and kinetic energy. Additionally, 
the good a greement between th e high-resolution low-7 12 CO ob- 
servations (IZapata et all 120091) and our higher-7 12 CO data indi- 
cates a close spatial correlation between the IR H2 fingers and 
the molecular outflows, which are likely driven by an explosive 
event instead of the outflows from a single young stellar object or 
radio source, i.e., SMA-1 or IRc2. Nevertheless, if this explosive 
outflow is a superposition of several flows emanating from pro- 
tostellar objects in BN/KL, they are likely distributed uniformly 
in position angle and will show little asymmetry in redshifted 
and blueshifted components, which are ruled out by our results. 



IZapata et al.1 (120091 120111) showed that the 12 CO bullets are 
less present in the southern and southeastern directions (see also 
Fig.[3]b). Because the N-S dense ridge in Orion BN/KL has a 
velocity gradient where the northern part (with respect to the 
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5 h 35 m 16 s 14 s 
1 R.A. (J2000) 

Fig. 4. SMA high-velocity 12 CO bullets overlaid with the JVLA 
methanol maser emission (K. M. Menten, priv. comm.). The red 
and blue 12 CO emission contours are the same as shown in Figs. 
[T]and [3] The CH3OH 42,2 - 4i 3 E maser emission is shown in 
grey contours plotted 0.05 to 3.05 Jy beam -1 in step of 0.2 Jy 
beam -1 . The black dashed line indicates the possible symmet- 
ric distribution between the I2 CO and CH3OH emission, where 
most of the high-velocity 12 CO bullets are likely absorbed below 
this line. The positions of the explosion center (EC), source BN, 
HC, and H2 peak 1 and 2 are marked. Magenta triangles repre- 
sent the positions of source /, SMA-1, source n, and IRc7 (east 
to west). The three stron gest methyl for mate emission peaks are 
marked in black stars dFavre et all 1201 ll) . 



HC) has a LSR veloc ity of > 9 km s" 1 , and < 8 km s -1 for the 
denser southern part dPeng et all 1201 2|) . it is likely that most of 
the bullets are decelerated and eventually absorbed in the south- 
ern part of the cloud. This scenario is supported by the detection 
of methanol masers (e.g., the 42,2 - 4i 3 E line shown in Fig.|4|, 
which are located mostly in the southern part o f the ridge and 
appea r to have formed behind the shock fronts ( Joh nston et all 
[1991 . Therefore, it is possible that the shocks due to this explo- 
sive outflow are absorbed in the southern part of Orion BN/KL 
and induce the methanol masers. There may be temperature dif- 
ferences in the cloud that affect the formation and/or distribution 
of some molecules. This scenario may help to explain the long- 
debated different iation between oxygen- and nitrogen-bearing 
species (see, e.gJBlake et aTlll987tlRodgers & Charnievl fizOOl; 
iLiu et all 120021 iFriedel & Snvder. l2008l) . namely that its main 
cause may be related to shocks. 

Furthermore, if one- forth of the bullets (~ 2.5 x 10 46 erg) 
are absorbed by the cloud in 500 years, the resulting luminosity 
(L ~ E/t) is > 400 L , where the time scale for absorption and 
re-radiation is shorter than 500 years because the ejecta would 
take time to reach the outer part of the cloud. T he total bolo- 
metri c luminosity for BN/KL is about 7 x 10 4 L dGezari et all 
1998) assuming a distance of 414 pc, and therefore the ex- 
plosive outflow may not be energetic enough to illuminate the 
whole BN/KL region (other stellar objects may also contribute). 
However, IRc2 with a luminosity of just 1000+5 00 L n and the 
other less luminous infrared sources in BN/KL dGezari et all 
1998) may be illuminated by this outflow. 



4. Summary 

We detected very high velocity (up to ±100 km s -1 ) outflows 
in the 12 CO 7=6-5 and 7=7-6 transitions in the Orion BN/KL 
region. The morphology of our 12 CO emission agrees well with 
the high-resolution SMA observations, where the molecular out- 
flow is driven by an explosive event and is related to the IR H2 
fingers seen on a larger scale in Orion BN/KL. The kinetic en- 
ergy of this explosive outflow calculated from the high-velocity 
CO bullets is estimated to be about 1 - 2 x 10 47 erg. It is still 
not clear what caused this explosion, but it is likely related to 
the runaway objects from a dynamically decayed multiple sys- 
tem in Orion BN/KL. We also proposed a scenario in which the 
methanol masers are induced by the same outflow at which most 
of the shocks are decelerated and absorbed in the southern part 
of BN/KL. 
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